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This is a continuation of two earlier papers on the subject of reflection of waves from 
inhomogeneous isotropic media. In this particular paper an exponential conductivity 
profile is perturbed in such a manner that the conductivity is increased for all heights above 
a certain level. A hyperbolic tangent transition is employed in order to avoid discontinuities 
in the conductivity versus height profile. 



1. Introduction 

In two earlier communications from the present 
authors [Wait and Walters 1963a and L 963b] oblique 
reflection of (VLF) radio waves from continuously 
stratified ionized media was considered. In (he 
first of these, denoted (I), the profile of (he effective 
conductivity was taken to be exponential in form. 
It was pointed out that an isotropic exponential 
model is a fairly good representation of the I) layer 
of the ionosphere under quiet DAYTIME condit ions. 
In the second paper, denoted (II), an idealized 
perturbed exponential model was considered. The 
perturbation consists of a localized increase of elec- 
tron density which itself has a Gaussian profile. The 
reflection coefficient was shown to he influenced by 
the vertical location of this Gaussian perturbation. 
In the present paper, which is part 111 of the series, 
the exponential profile model introduced in (I), is 
modified by allowing the electron density to he 
increased for all height s above a specified level. Such 
a modification to the quiescent ionosphere could 
result from ionizing radiations associated with a solar 
flare [Jean and Crary, 1962; Pierce, 1963; Chilton 
et aL, 1963]. 

For the detailed numerical results given here, it is 
assumed that the earth's magnetic field may be 
neglected. This is justified for highly oblique 
incidence at VLF provided that attention is restricted 
to effects which result from ionization in lowest 
da}- time ionosphere [J older, 1962]. 

2. Description of the Profile 

The notation follows that used in (I) as closely 
as possible. Thus, the undisturbed profile, as a 
function of height z, is defined by the conductivity 
parameter 1/L(z) where 

1/Z(z)=(l/Z) exp (0z), (1) 

and L is a constant, (3 is a gradient parameter, and 



z is the height above the reference level. Under 

the isotropic assumption, it is known that [Wait, 
1962] 

Z= ^M, (2 ) 



in terms of the angular frequency co, collision fre- 
quency v, and the plasma frequency co . Further- 



more, at VLF, y>>co which leads to 
L^o)/u r where co r =ul/v. 



(3) 



In general, the conductivity parameter is propor- 
tional to N(z)/v(z) where N(z) and v(z) are the 
electron density and collision frequency regarded as 
a function of height. The constant 13, in the ex- 
ponent, is a measure of the sharpness of the gradient. 
As in (II), it will be assumed that (8=0.3 km -1 
typifies an undisturbed /Mayer profile. 

We shall now turn our attention to the modifica- 
tion of the exponential profile. As in (II), it is 
assumed that the collision frequency profile is 
unchanged whereas the ionization is to be increased 
by an amount AN(z) where 

AiV(s)=AA^ [tanh(^)+l} (4) 

and ATVo, F, and D are const aids. It is evident 
that AN(z) is equal to AN n at z=F which may be 
described as the "lower edge" of the modified layer. 
When z is somewhat less than F it is seen that 
AN(z) is zero whereas when z is much greater than 
F, AN(z) becomes equal to 2A7V . It is clear that 
the vertical location of the "lower edge" is governed 
by the value of F. The rapidity of the transition 
from to 2AA r o is controlled by the magnitude of D; 
the quantity 21) could be described as the "transi- 
tion thickness." 

Following the discussion given in (II), the collision 
frequency profile is taken to be 



v(z)=v Q exp (— j8s/2), 



(5) 
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where 0=0.30 km \ The resulting conductivity 
perturbation has the form 



AN(z) ANq 



Vq 



exp Gfe/2) [l + tanh (^)} (6) 



The complete conductivity profile is now given by 



L{z) L(0) 



/exp03 2 )+.4exp(|s) 



[> 



-tanh ( 



z-F 
V D 



)]}• 



(7) 



where the right-hand side is proportional to the 
effective conductivity of the medium as a function 
of height above (or below) the reference level at 
z=0. The coefficient A defines the strength or 
magnitude of the perturbation. The constant i(0) 
is chosen to be equal to 7.5/X where X is the wave- 
length in kilometers. 

To enable the reader to grasp the significance of 
the various possible forms of profiles, some typical 
examples are illustrated in figures la, lb, and lc. 
In figure la, the vertical location of the "lower edge" 
of the perturbed layer is varied while other factors 
are kept constant. In lb, the abruptness of the 
transition at the "lower edge" is varied for fixed 
values of F and A. Finally, in lc, the strength of 
the perturbed layer is varied while keeping F and D 
constant. 

3. Results of the Calculations 

The method used in the reflection coefficient cal- 
culations is the same as that used in (I) and (II), 
so this aspect of the subject need not be discussed. 
As before, the amplitude \R\ and the phase of R 
for a vertically polarized wave are considered. 
The reflection coefficient is evaluated in the free- 
space region corresponding to 2^— °o; however, 
it is to be remembered that the phase is referred to 
the level 2=0. 

In figure 2a, the amplitude of the reflection coeffi- 
cient is plotted as a function of F for X=15 km 
(i.e., /=20 kc/s), A=2, D=2 km, and C values 
varying from 0.05 to 0.4. Small values of C cor- 
respond to angles near grazing; these are important 
for long distance propagation of VLF radio waves. 
It is evident that when F is somewhat above zero, 
\R\ assumes the value for an unperturbed exponential 
layer. As the "lower edge" of the perturbed layer 
moves down (i.e., F becoming negative), the tend- 
ency is for \R\ to increase somewhat. As the per- 
turbed layer is moved farther down (i.e., becoming 
increasingly negative), \R\ tends to be reduced. 
The behavior illustrated here is not dissimilar to 
the corresponding case of a Gaussian "bump," 
described in (II), as it moves down to low heights. 
In both cases, the perturbation enhances the re- 
flection at greater heights but degrades it at lower 
heights. 

The phase of R is shown in figure 2b for the same 
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Figures 2a, 2b. Amplitude and phase of the reflection coefficient 
for various values of C, the cosine of the angle of incidence. 
F is the vertical distance, relative to the reference height, to the "lower edge" 
of the perturbed layer.] 
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Figures 3a, 3b. Amplitude and phase of the reflection coefficient 
for various values of D, the "half .thickness" of the layer 
transition for C = 0.2. 
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Figures 4a, 4b. Amplitude and phase of the reflection coefficient 
for various values of D, the "half -thickness" of the layer 
transition for G = 0. 1 . 
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Figures 5a, 5b. Amplitude and phase of the reflection coefficient 
for various values of X, the wavelength. 
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Figures 6a, 6b. Amplitude and phase of the reflection coefficient 
for various values of A. the magnitude of the perturbation. 



conditions as in figure 2a. It is apparent that when 
the "lower edge" is somewhat above F=Q, the phase 
attains its undisturbed value. As F becomes neg- 
ative the phase increases and, for grazing angles, 
this increase is mono tonic in nature. However, for 
the steeper angles of incidence there is some evidence 
of oscillations, which are presumably due to certain 
interference phenomena. 

The dependence of the "transition thickness" 
on the reflection coefficient is shown in figures 3a 
and 3b for (7=0.2. These show, as one might ex- 
pect, that a sharp transition is associated usually 
with an increased reflection coefficient. As the 
transition is stretched out, the tendency is for the 
maximum at 2^ — 15 km to disappear. It is 
particularly interesting to note that the phase is not 
overly sensitive to changes of D. 

The influence of varying I) for a highly grazing 
situation is shown in figures 4a and 4b. These 
show behavior which is qualitatively the same as 
figures 3a and 3b. 

The wavelength dependence of the reflection 
coefficient is shown in figures 5a and 5b for A = 10, 
15, 20, 25, and 30 km corresponding to frequencies of 
30, 20, 15, 12, and 10 kc/s. As expected, the curves 
have similar shapes to each other and the tendency 
is for the shorter wavelengths to be associated with 
more pronounced changes. 

Finally, in figures 6a and 6b, the influence of 
varying the strength A of the perturbed layer is 
indicated. Generally, as A is increased, the changes 
are more pronounced. It is interesting to note that 
at the very low heights the magnitude of the re- 
flection coefficient is very sensitive to small changes 
of A. This behavior is similar to that found in the 
troposphere when considering the reflection of VHF 
radio waves from gradient changes of the refractive 
index [Wait, 1962], 

4. Concluding Remarks 

In the present series of three papers, the reflection 
of plane waves from an inhomogeneous isotropic 
medium has been considered. Attention has been 
confined to horizontally stratified layers whose 
effective conductivity varies smoothly with height. 
Special attention is given to vertically polarized 
waves at highly oblique incidence and at very low 
frequencies because these are important in long 
distance propagation. 

In the first paper, the conductivity profile is 
taken to be exponential in form. Although this is 
a simple model, it is quite convenient to employ for 
quantitative estimates of the amplitude and phase 
changes as a function of gradient, angle of incidence, 
and frequency. The application of the results of 
(I) to mode theory has already been carried out 
[Wait, 19631. 

In the second paper of this series, denoted (II), 
the exponential profile is allowed to have a perturba- 
tion in the form of a Gaussian curve. In the third 
paper (i.e., the present one), the perturbation is in 
the form of a hyperbolic tangent. Basically, the 
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difference between these two classes of profiles is 
that in (II) the perturbation is localized whereas in 
(III) it is spread out over a range of heights. It is 
probable that the results in (II) are more applicable 
to ionization produced by the hard X-rays emitted 
from a space nuclear burst. On the other hand, the 
results in (III) would be more appropriate to the 
ionization produced by a solar flare. 

The direct application of the results in (II) and 
(III) to mode theory computations and further 
analysis including the earth's magnetic field is 
reserved for future communications. In the mean- 
time, it is hoped that the present results should add 
to the knowledge of wave propagation in inhomo- 
geneous media. 

We thank A. G. Jean for many useful suggestions 
concerning this work. The work has been supported 
by the Advanced Research Projects Agency, Wash- 
ington, D.C., under ARPA Order No. 183-62. 
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